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X-RAY, OPTICAL, AND RADIO OBSERVATIONS OF THE TYPE II SUPERNOVAE 1999EM AND 1998S 
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ABSTRACT 

Observations of the Type II-P (plateau) Supernova (SN) 1999em and Type Iln (narrow emission line) SN 1998S 
have enabled estimation of the profile of the SN ejecta, the structure of the circumstellar medium (CSM) estab- 
lished by the pre-SN stellar wind, and the nature of the shock interaction. SN 1999em is the first Type II-P 
detected at both X-ray and radio wavelengths. It is the least radio luminous and one of the least X-ray luminous 
SNe ever detected (except for the unusual and very close SN 1987A). The Chandra X-ray data indicate non- 
radiative interaction of SN ejecta with a power-law density profile {p oc r~" with n ^1) for a pre-SN wind with 
a low mass-loss rate of ^ 2 x 10"'' Mq yr"' for a wind velocity of 10 km s"', in agreement with radio mass-loss 
rate estimates. The Chandra data show an unexpected, temporary rise in the 0.4-2.0 keV X-ray flux at ~100 days 
after explosion. SN 1998S, at an age of >3 years, is still bright in X-rays and is increasing in flux density at cm 
radio wavelengths. Spectral fits to the Chandra data show that many heavy elements (Ne, Al, Si, S, Ar, and Fe) 
are overabundant with respect to solar values. We compare the observed elemental abundances and abundance 
ratios to theoretical calculations and find that our data are consistent with a progenitor mass of approximately 
15-20 Mq if the heavy element ejecta are radially mixed out to a high velocity. If the X-ray emission is from 
the reverse shock wave region, the supernova density profile must be moderately flat at a velocity ^ lO'* km s"', 
the shock front is non-radiative at the time of the observations, and the mass-loss rate is 1-2x10"'* yr"' for 
a pre-supernova wind velocity of 10 km s"'. This result is also supported by modeling of the radio emission 
which implies that SN 1998S is surrounded by a clumpy or filamentary CSM established by a high mass-loss rate, 
' 2 X 10"^ M(7i yr"', from the pre-supernova star. 
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(SN 1998S, SN 1999em) — X-rays: ISM — radio continuum: ISM 



X-rays: individual 



Suzuki & Nomoto 1995; Zimmermann et al. 1994), 19941 
(Immler, Pietsch, & Aschenbach 1998b; Immler, Wilson, & 
Terashima 2002), 1994W (Schlegel 1999b), 1995N (Fox et al. 
2000), 1998bw (aka GRB 980425; Plan et al. 1999; Plan et al. 
2000), 1999em (Fox et al. 1999; Schlegel 2001a), and 1999gi 
(Schlegel 2001b). We present here results from Chandra ob- 
servations of SN 1999em and SN 1998S. In general, the high 
X-ray luminosities of these 14 SNe (Lx lO^^-lO"*' ergs s"') 
dominate the total radiative output of the SNe starting at an age 
of about one year. This soft X-ray (< 10 keV) emission is most 
convincingly explained as thermal radiation from the "reverse 
shock" region that forms within the expanding SN ejecta as it 
interacts with the dense stellar wind of the progenitor star 
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1. INTRODUCTION 

To date, 14 supernovae (SNe) have been detected in X-rays 
in the near aftermath (days to years) of their explosions'"*: 
SNe 1978K (Peti-e et al. 1994; Schlegel 1995; Schlegel, Pe- 
tre, & Colbert 1996; Schlegel 1999a), 1979C (Immler, Pietsch, 
& Aschenbach 1998a; Kaaret 2001; Ray, Petre, & Schlegel 
2001), 1980K (Canizares, Kriss, & Feigelson 1982; Schlegel 
1994, 1995), 1986J (Houcket al. 1998; Schlegel 1995; Chugai 
1993; Bregman & Pildis 1992), 1987 A (Burrows et al. 2000; 
Dennerl et al. 2001; Gorenstein, Hughes, & Tucker 1994; 
Hasinger, Aschenbach, & Triiumper 1996), 1988Z (Fabian & 
Terlevich 1996), 1993J (Immler, Aschenbach, & Wang 2001; 
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The interaction of a spherically symmetric SN shock and a 
smooth circumstellar medium (CSM) has been calculated in de- 
tail (Chevalier & Fransson 1994; Suzuki & Nomoto 1995). As 
the SN "outgoing" shock emerges from the star, its characteris- 
tic velocity is ^lO'* km s"', and the density distribution in the 
outer parts of the ejecta can be approximated by a power- law 
in radius, p oc r~", with 7 < n < 20. The outgoing shock prop- 
agates into a dense CSM formed by the pre-SN stellar wind. 
This wind is slow (v„ ~ 10 km s"') and was formed by a pre- 
supernova mass-loss rate of M ~ 10""* to 10"^ yr"'. The 
CSM density for such a wind decreases as the square of the ra- 
dius {p = M/47rr^Vw). The colUsion between the stellar ejecta 
and the CSM produces a "reverse" shock, which travels out- 
ward at '^lO^ km s"' slower than the fastest ejecta. Interaction 
between the outgoing shock and the CSM produces a very hot 
shell (~10' K) while the reverse shock/ejecta interaction pro- 
duces a denser, cooler shell (~10^ K) with much higher emis- 
sion measure and is where most of the observable X-ray emis- 
sion arises. If either the CSM density or n is high, the reverse 
shock is radiative, resulting in a dense, partly-absorbing shell 
between the two shocks. 

Chevalier (1982) proposed that the outgoing shock from the 
SN explosion generates the relativistic electrons and enhanced 
magnetic field necessary for synchrotron radio emission. The 
ionized CSM initially absorbs most of this emission, except 
in cases where synchrotron self-absorption dominates, as may 
have been the case in SN 1993J (Chevalier 1998; Fransson & 
Bjornsson 1998). However, as the shock passes rapidly out- 
ward through the CSM, progressively less ionized material is 
left between the shock and the observer, and absorption de- 
creases rapidly. The observed radio flux density rises accord- 
ingly. At the same time, emission from the shock region is 
decreasing slowly as the shock expands so that, when radio ab- 
sorption has become negligible, the radio light curve follows 
this decline. Observational evidence also exists at optical wave- 
lengths for this interaction of SN shocks with the winds from 
pre-supernova mass loss (Filippenko 1997). 

All known radio supernovae appear to share common proper- 
ties of (i) non-thermal synchrotron emission with a high bright- 
ness temperature, (ii) a decrease in absorption with time, re- 
sulting in a smooth, rapid turn-on first at shorter wavelengths 
and later at longer wavelengths, (ui) a power-law dechne of the 
emission flux density with time at each wavelength after max- 
imum flux density (absorption r « 1) is reached at that wave- 
length, and, (iv) a final, asymptotic approach of the spectral in- 
dex a to an optically thin, non-thermal, constant negative value 
(Weiler et al. 1986). These properties are consistent with the 
ChevaUer model. 

The signatures of circumstellar interaction in the radio, opti- 
cal, and X-ray regimes have been found for a number of Type 
II SNe such as die Type II-L SN 1979C (Fesen & Maton- 
ick 1993; Inmiler, Pietsch, & Aschenbach 1998a; Weiler et al. 
1986, 1991) and tiie Type II-L SN 1980K (Canizares, Kriss, & 
Feigelson 1982; Fesen & Becker 1990; Leibundgutet al. 1991; 
Weiler et al. 1986, 1992). The Type lln subclass has peculiar 
optical characteristics: narrow Ha emission superposed on a 
broad base; lack of P Cygni absorption-hne profiles; a strong 
blue continuum; and slow evolution (Schlegel 1990; Fihppenko 
1997). The narrow optical lines are clear evidence for dense 
circumstellar gas — they probably arise from the reprocessing 
of the X-ray emission — and are another significant means by 
which the shock radiatively cools. The best recent examples 



of Type lln SNe include SNe 1978K, 1986J, 1988Z, 1994W 
(which may have been a peculiar Type 11-P, Sollerman, Cum- 
ming, & Lundqvist 1998), 1995N, 1997eg, and 1998S. 

SN 1999em was opticaUy discovered on 1999 October 29 in 
NGC 1637, approximately 15" west and 17" south of the galac- 
tic nucleus (Li et al. 1999). The spectrum on 1999 October 29 
shows a high temperature (Baron et al. 2000), indicating that 
the supernova was discovered at an early phase. We take 1999 
October 28 as the date of the explosion. SN 1999em reached 
a peak brightness of my « 13.8 mag on about day 4, and re- 
mained in the "plateau phase," an enduring period of nearly 
constant V-band brightness, until about 95 days after discovery 
(Leonard et al. 2001b). At an estimated distance of 7.8 Mpc 
(Sohn & Davidge 1998; see also Hamuy et al. 2001 who find 
a distance of 7.5 Mpc to SN 1999em based on the Expanding 
Photosphere Method), this is the closest Type II-P supernova 
yet observed, and it is the only Type II-P observed in both X- 
rays and radio. Its maximum brightness of My ~ -15.75 marks 
it as one of the optically least luminous SNe 11-P (cf. Patat 
et al. 1994). Smartt et al. (2001) recently derived a mass of 
12 ± 1 Mq for file progenitor of SN 1999em. 

On 1998 March 3, SN 1998S was optically discovered 16" 
west and 46" south of the nucleus of NGC 3877 (Li et al. 1998). 
Its spectrum quickly revealed it to be a peculiar Type Iln, with 
narrow H and He emission lines superimposed on a blue con- 
tinuum (Filippenko & Moran 1998) along with strong N III, 
C II, C III, and C IV emission (Garnavich, Jha, & Kirshner 
1998a; Liu et al. 2000) similar to the "Wolf-Rayet" features 
seen in the early observations of SN 1983K (Niemela, Ruiz, & 
Phillips 1985). SN 1998Sreachedapeak brightness of my » 12 
(Granslo et al. 1998; Garnavich et al. 1998b). For a distance of 
17 Mpc (Tully 1988), the intrinsic brightness of My « -18.8 
marks SN 1998S as one of the optically brightest Type II SNe 
(cf. Patat et al. 1993). 

2. OBSERVATIONS 
2.1. SN1999em 

2.1.1. X-ray 

The trigger for our first Chandra observation was the optical 
discovery of a SN within 10 Mpc. The extremely rapid response 
of the Chandra staff allowed us to catch SN 1999em three days 
after discovery. We observed SN 1999em on 1999 November 
1, November 13, December 16, 2000 February 7, October 30, 
2001 March 9, and July 22 (days 4, 16, 49, 102, 368, 495, and 
630 from reference). A summary of the observations is given in 
Table 1, and Schlegel (2001a) also discusses the first three. All 
observations were performed by the Advanced CCD Imaging 
Spectrometer (ACIS) with the telescope aimpoint on the back- 
side illuminated S3 chip, which offers increased sensitivity to 
low energy X-rays over the front-side illuminated chips. The 
first two observations used a 1/2 subarray mode (only half of 
the CCD was read out), while the rest used the full CCD. The 
data were taken in timed-exposure (TE) mode using the stan- 
dard integration time of 3.2 sec per frame and telemetered to 
the ground in very faint (VF) mode for the first two observa- 
tions and faint (F) mode for the rest. In VF mode, the teleme- 
tered data contain the values of 5 x 5 pixel islands centered on 
each event, while in F mode the islands are 3 x 3 pixels. 

2.1.2. Radio 
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We observed SN 1999em over 30 times with the Very Large 
Array (VLA)'^ at frequencies ranging from 43.3 GHz to 1.5 
GHz and ages from 3 to 454 days after the discovery date of 
1999 October 29 (Lacey et al. 2002). Except for six weak but 
significant (> 3(t) detections at mid-cm wavelengths (1.5 - 8.4 
GHz) at ages of between 33 and 69 days, all measurements 
yielded only upper Umits. 

2.1.3. Optical 

Extensive optical ground-based observations were made of 
SN 1999em; these include f/BVT?/ photometry, long-term spec- 
troscopy, and spectropolarimetry (Leonard et al. 2001a). These 
data were used to obtain a distance with the Expanding Pho- 
tosphere Method (Leonard et al. 2001b). Both the photometry 
and spectroscopy suggest that this was a relatively normal Type 
11-P event. We note, however, that SN 1999em was about 1 
magnitude fainter in the V-band during the plateau than the av- 
erage of the 8 SNe II-P with published photometry and previ- 
ously derived EPM distances, and it may have also had a some- 
what unusual color evolution (Leonard et al. 2001b). 

Spectropolarimetry of SN 1999em was obtained on 1999 
November 5, December 8, and December 17 (7, 40, and 49 
days after discovery, respectively, while it was still in its plateau 
phase) with the Kast double spectrograph (Miller & Stone 
1993) with polarimeter at the Cassegrain focus of the Shane 
3-m telescope at Lick Observatory. Similarly, on 2000 April 5 
and 9 (159 and 163 days after discovery, long after the plateau). 
The Low-Resolution Imaging Spectrometer (Oke et al. 1995) 
was used in polarimetry mode (LRISp; Cohen 1996) at the 
Cassegrain focus of the Keck-1 10-m telescope. Leonard et al. 
(2001a) discuss the details of the polarimetric observations and 
reductions. 

Total flux spectra were obtained at many epochs, primar- 
ily with the Nickel 1-m and Shane 3-m reflectors at Lick Ob- 
servatory. UBVRI photometry was conducted with the 0.8-m 
Katzman Automatic Imaging Telescope (KAIT; Li et al. 2000; 
Filippenko et al. 2001) at Lick Observatory. These data and 
their implications are thoroughly discussed by Leonard et al. 
(2001b). 

2.2. SN1998S 

2.2.1. X-ray 

Our criterion for Chandra observations of SNe beyond 
10 Mpc was the detection of radio emission greater than 1 mJy 
at 6 cm. SN 1998S met this criterion on 1999 October 28 (Van 
Dyk et al. 1999), and subsequent reanalysis of early radio ob- 
servations indicated very weak detection at 8.46 GHz as early as 
1999 January 07 (Lacey et al. 2002). We observed with Chan- 
dra on 2000 January 10, March 7, August 1, 2001 January 14, 
and October 17 (days 678, 735, 882, 1048, and 1324 since opti- 
cal discovery). Similar to the observations of SN 1999em, these 
data were taken with the AC1S-S3 chip at nominal frame time 
(3.2 sec) in TE mode and telemetered in F mode. A summary 
of the observations is Usted in Table 2. 

2.2.2. Radio 

We observed SN 1998S over 20 times with the VLA at fre- 
quencies ranging from 22.48 GHz to 1.46 GHz and ages from 

The VLA is operated by the NRAO of the AUI under a cooperative agreement with the NSF. 

These spectra are produced from all available source photons, as opposed to the "polarized flux spectra" which refer only to the net polarization as a function of 
wavelength. 



89 to 1059 days after the discovery date of 1998 March 03 
(Lacey et al. 2002). All observations yielded only upper lim- 
its until 1999 January 07, age 310 days, when a weak detection 
was obtained. Since that time, monitoring has continued and 
SN 1998S has increased in flux density at all mid-cm wave- 
lengths. 

2.2.3. Optical 

On 1998 March 7, just 5 days after the discovery of 
SN 1998S, optical spectropolarimetry was obtained of the su- 
pernova with LRISp on the Keck-II 10-m telescope (Leonard 
et al. 2001a). Total flux spectra were obtained over the first 
494 days after discovery, as follows: (a) at Keck using LRIS on 
1998 March 5, 6, and 27, and on 1999 January 6, and (b) at flie 
Lick 3-m Shane reflector using the Kast double spectrograph on 
1998 June 18, July 17 and 23, and November 25, and on 1999 
January 10, March 12, and July 7. See Leonard et al. (2000) for 
details of the optical observations and reductions. 

3. Chandra data reduction 

For each observation, we followed the data preparation 
threads provided by the Chandra X-ray Center. We used the 
Chandra Interactive Analysis of Observations (ClAO) soft- 
ware (version 2.2) to perform the reductions, along with the 
CALDB 2.10 calibration files (gain maps, quantum efficiency, 
quantum efficiency uniformity, effective area). Bad pixels were 
excluded, and intervals of bad aspect were eliminated. 

All data were searched for intervals of background flaring, 
in which the count rate can increase by factors of up to 100 
over the quiescent rate. The Ught curves of background regions 
were manually inspected to determine intervals in which the 
quiescent rate could be reliably calculated. In accordance with 
the prescription given by the ACIS Calibration Team, we iden- 
tify flares as having a count rate >30% of the quiescent rate. 
Background flares were found in each of the first four and the 
last observation of SN 1999em, lasting from hundreds of sec- 
onds to over 15 ksec. Only the second and fifth observations 
of SN 1998S had such flares, lasting ~300 sec and --14 ksec, 
respectively. The event lists were filtered to exclude the time 
intervals during which flares occurred. Effective observation 
times after filtering are listed in Tables 1 and 2. 

In addition, we have performed a reprocessing of all data so 
as not to include the pixel randomization that is added during 
standard processing. This randomization has the effect of re- 
moving the artificial substructure (Moire pattern) that results as 
a byproduct of spacecraft dither Since all of our observations 
contained a substantial number of dither cycles (one dither cy- 
cle has a period of ~1000 sec), this substructure is effectively 
washed out, and there is no need to blur the image with pixel 
randomization. Removing this randomization slightly improves 
the point spread function (PSF). 

The source spectra were extracted with the CIAO tool dmex- 
tract. We excluded events with a pulse invariant (PI) of either 
(underflow bins) or 1024 (overflow bins). The region of extrac- 
tion was determined from the CIAO tool wavdetect, a wavelet- 
based source detection algorithm which can characterize source 
location and extent. In each observation, the region determined 
by wavdetect for the SN is consistent with that of a point source. 
The data were restricted to the energy range of 0.4-8.0 keV 
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Fig. 1. — NGC 1637 in optical and X-rays. X-ray contours from the first Chandra observation of SN 1999em are overlaid on a DSS image of the host galaxy. 



for the purposes of spectral analysis since the effective area of 
ACIS falls off considerably below 0.4 keV and the increased 
background above 8 keV makes this data unreliable. 



observation to the fourth, despite the continued decline of the 
high-energy X-rays. The multicolour X-ray Ughtcurve is shown 
in Fig. 2. 



4. RESULTS 

4.1. SN 1999em 
4.1.1. X-my 

The source was detected on all occasions, and an overlay 
of the X-ray contours from the first observation on an opti- 
cal image of the host galaxy is shown in Fig. 1. The wavde- 
tect position" (based solely on the Chandra aspect solution) 
is 4Hl"'27!'16, -2°51'45."6, which is within 1" (roughly the 
Chandra pointing accuracy) of the radio and optical positions. 
The low flux (/x < 10"^* ergs cm~^ s~^) resulted in few to- 
tal counts (see Table 1). We have fit both the MEKAL and 
thermal Bremsstrahlung models to our data. The MEKAL 
model is a single-temperature, hot diffuse gas (Mewe, Gro- 
nenschild, & van den Oord 1985) with elemental abundances 
set to the solar values of Anders & Grevesse (1989). We have 
calculated the column density (Predehl & Schmitt 1995) to be 
A'h = 6. 1 X 10^" cm"^ based upon a value of E(b-v) = 0.05 (Baron 
et al. 2000). For both models, only the temperatures and over- 
all normalizations were allowed to vary. In XSPEC (Arnaud 
1996), we performed maximum likelihood fits on the unbinned 
data using Cash statistics (Cash 1979). These fits should be 
insensitive to the number of counts per channel and are thus ap- 
propriate for fitting low-count data. The best-fit observed tem- 
peratures are shown in Table 3. As expected, the uncertainties 
in fit parameters for such low-count data are large. 

We have also used flux ratios to characterize the source. The 
fluxes were calculated by manually integrating the source spec- 
tra, i.e., they are model-independent, absorbed fluxes. The 
first observation had a flux ratio of /x[2-8 keV]//x[0.4-2 keV] 
= 2.1 ±0.9, but later observations were softer (Table 1). Some- 
what surprisingly, the total flux nearly doubled from the third 

This is the average position of observations 2-6, in which the standard deviations 
was off by ~3". This offset was corrected for Fig. 1. 



4.1.2. Radio 

After 15 unsuccessful VLA observations of SN 1999em at 
ages 4 to 19 days after the assumed explosion date of 1999 Oc- 
tober 28, the SN was finally detected on day 34 with 0.190 mJy 
flux density at 8.435 GHz. Subsequent observing on days 46, 
60, and 70 indicated that SN 1999em remained detectable at 
a similar flux density, near the sensitivity limit of the VLA, at 
8.435, 4.885, and 1.465 GHz. The data are shown, together 
with the x-ray results, in Fig. 2. Since day 70, no radio detection 
of SN 1999em in any cm wavelength band has been obtained. 
SN 1999em is among only a few radio detected Type II-P, and 
its behavior is difficult to compare with other radio observa- 
tions of the same type SN. An estimated spectral luminosity of 
SN 1999em at 6 cm peak on day ~34 of Lgcm peak ~ 2.2 x 10^^ 
erg s"' Hz ' makes it the least radio luminous RSN known 
except for the peculiar, and very nearby, SN 1987 A. The ra- 
dio position of the emission from SN 1999em is 4*'41™27.^157, 
-2°51'45."83 (J2000). The details of the observations and anal- 
ysis will be found in Lacey et al. (2002). 

4.1.3. Optical 

The spectropolarimetry of SN 1999em provided a rare op- 
portunity to study the geometry of the electron-scattering at- 
mosphere of a "normal" core-coUapse event at multiple epochs 
(Leonard et al. 2001a). A very low but temporally increasing 
polarization level suggests a substantially spherical geometry at 
early times that becomes more aspherical at late times as ever- 
deeper layers of the ejecta are revealed. When modeled in terms 
of oblate, electron-scattering atmospheres, the observed polar- 
ization implies an asphericity of at least 7% during the period 
studied. We speculate that the thick hydrogen envelope intact 

were ctra = 0"03 and croec = 0"4. Observation 1 had a known aspect error and 
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Fig. 2. — X-ray and radio light curves of SN 1999em. The detections are shown with error bars while the upper limits are shown with downward arrows. Radio 
flux densities are plotted with the appropriate symbols. The absorbed X-ray flux is plotted in both the 2-8 keV (solid line) and 0.4-2 keV (dashed line) bands. The 
hard X-rays follow a steady decline, but the soft X-rays double from the third observation to the fourth. 



at the time of explosion in SNe II-P might serve to dampen the 
effects of an intrinsically aspherical explosion. The increase 
in asphericity seen at later times is consistent with a trend re- 
cently identified (Wang et al. 2001) among stripped-envelope 
core-collapse SNe: the deeper we peer, the more evidence we 
find for asphericity. The natural conclusion that it is an explo- 
sion asymmetry that is responsible for the polarization has fu- 
eled the idea that some core-collapse SNe produce gamma-ray 
bursts (GRB; e.g., Bloom et al. 1999) through the action of a 
jet of material aimed fortuitously at the observer, the result of 
a "bipolar," jet-induced, SN explosion (Khokhlov et al. 1999; 
Wheeler et al. 2000; Maeda et al. 2000). Additional multi- 
epoch spectropolarimetry of SNe 11-P are clearly needed, how- 
ever, to determine if the temporal polarization increase seen in 
SN 1999em is generic to this class. 

4.2. SN1998S 
4.2.1. X-ray 

The source was detected on all occasions, and an overlay 
of the X-ray contours from the first observation on an opti- 
cal image of the host galaxy is shown in Fig. 3. The wavde- 
tect position^^ (based solely on the Chandra aspect solution) 
is llHe^OefM, 47°28'55."1 (J2000), which is within 0"5 of 
the radio and optical positions. In addition to constructing 
multicolour X-ray lightcurves (Fig. 4), the high luminosity of 
SN 1998S (Lx « 10^0 ergs s"') aUowed for basic spectral fitting 
to be done for each of the first four observations; in the fifth, 
there were too few counts for a similar type of analysis. We 
performed this in XSPEC with the VMEKAL model, which 
is identical to the MEKAL model but allows individual ele- 
mental abundances to vary. Using the reddening E(b^v) = 0.23 
(Leonard et al. 2000), we calculated the hydrogen colunm den- 
sity to be A^H = 1 -36 X lO^' cm'^ (Predehl & Schmitt 1995). We 
used a redshift of z = 0.003 (Willick et al. 1997). 

Although there were sufficient counts (>500 per observa- 
tion) in each of the first four observations for spectral model- 
ing, there were few counts per energy channel. To mitigate the 
problems associated with fitting such data, we grouped the data 
to contain a minimum number of counts per bin. Because there 
were virtually no counts above ~ 7 keV in any of the observa- 
tions, we would "smear" the Fe line at 6.7 keV if we required 
more than ~ 7 counts per bin. Therefore, we set the minimum 

This is the average position of all observations. The standard deviations about 



number of counts per bin to the highest possible value which 
would still retain the integrity of the Fe line. This turned out to 
be 5 counts per bin for the first observation, 7 for the second, 5 
for the third, and 7 for the fourth. 

To reduce the number of free parameters in our fits, we froze 
the abundance of He to its solar value. We also linked the abun- 
dances of N, Na, Al, Ar, and Ca to vary with C. The remain- 
ing elemental abundances (O, Ne, Mg, Si, S, Fe, and Ni) were 
allowed to vary independently. We first allowed the column 
density to vary in our fits. Since the best-fit results were all 
consistent with the A^h obtained from the optical reddening, we 
fixed this parameter at the calculated value. 

We then performed ^ fitting using Gehrels weighting for the 
statistical error (Gehrels 1986), which is appropriate for low 
count data. We also performed maximum hkehhood fits to the 
unbinned data using Cash statistics (Cash 1979). The maxi- 
mum likelihood fits were in good agreement with the fits, 
and we report the results in Table 4. We note, however, that 
fitting with so few counts per bin may give misleading val- 
ues. For ease of viewing, we have plotted the best fit models 
against evenly binned data (Fig. 5). 

In order to bring out the line features, we have summed the 
individual spectra. We feel justified in doing this because the 
best-fit model parameters from the four individual observations 
with enough counts for individual modeling are consistent with 
each other, indicating that the spectral shape is not changing 
much, only the overall normalization (luminosity). We grouped 
the surmned spectrum to > 15 counts per bin for spectral fitting. 
Grouping to more counts per bin than that would have washed 
out some of the weaker line features. The individual response 
functions for each observations were combined in a weighted 
average based on exposure times. 

We fit four models to this summed spectrum: a single 
MEKAL plasma with solar abundances (1-T MEKAL); two 
MEKAL plasmas with solar abundances (2-T MEKAL); a 
single MEKAL plasma with variable abundances (all ele- 
ments allowed to vary except He; 1-T VMEKAL); and a sin- 
gle MEKAL with solar abundances plus a single MEKAL 
with variable abundances (2-T VMEKAL). The 1-T MEKAL 
model did not fit the data well, with a of 181 for 124 degrees 
of freedom (dof). Tests of the other three model fits against this 
one using the F-statistic indicate that only the 1-T VMEKAL 

means were (tra = 0"03 and ctdcc = 0" 17. 
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Fig. 3. — NGC 3877 in optical and X-rays. X-ray contours from the first Chandra observation of SN 1998S are overlaid on a DSS image of the host galaxy. 
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fit was a significant improvement (required at the 3.8 cr level). 
Most of the abundances in this model are not well constrained 
(especially those of C, N, Na, and Ni). The uncertainty in each 
abundance was determined in the following way. The abun- 
dance was changed from the best-fit value, the model was refit 
(with all other parameters allowed to vary), and the resultant 
was compared to the best- fit x^. This was repeated until 
the difference in was 2.706, which gave the 90% confidence 
interval for the given parameter. Table 5 lists the results. The 
best- fit temperature was kT = 9.8^j § keV, and the fit statistic 
was x^ = 139 for 116 degrees of freedom. Fig. 6 shows the 
summed spectrum and best-fit VMEKAL model with labels 
indicating the rough positions of elemental Unes. 

4.2.2. Radio 

Since its initial radio detection at 8.46 GHz on 1999 January 

07 (310 days after optical discovery), SN 1998S has increased 
at all mid-cm wavelengths. Modeling indicates that it should 
reach its 6 cm flux density peak at age ~1000 days, at roughly 
the epoch of the most recent observations on 2001 January 25. 
Assuming that is the case, SN 1998S is behaving like other 



radio-detected Type Hn (SNe 1986J and 1988Z both peaked at 
6 cm wavelength at age ~1000 days) although it will be the 
least radio luminous at 6 cm peak (Lfecm peak ^ 3.6 x 10^*" erg 
s-i Hz-i) of the six Type Iln RSNe known (SNe 1978K, 1986J, 
1988Z, 1995N, 1997eg, and 1998S). The radio data are shown 
together with the X-ray results in Fig. 4. 

Under the Chevalier model of a dense CSM established by a 
slow pre-supernova stellar wind (Vw = 10 km s"', Vshock = 10^ 
km s"\ r = 2 X 10"^ K; see, e.g., Weiler et al. 1986), such a ra- 
dio luminosity is interpreted as indicative of a mass-loss rate of 
^ 2x 10""* M0 yr~'. The radio position of the emission from 
SN 1998Sis ll''46'"6.n40,47°28'55."45. The details of the ob- 
servations and analysis will be found in Lacey et al. (2002). It 
should be noted that the best fit to the radio data requires signif- 
icant clumping or filamentation in the CSM, as was also found 
for the Type lln SNe 1986J and 1988Z (Weiler et al. 1990; Van 
Dyketal. 1993). 

4.2.3. Optical 

SN 1998S exhibits a high degree of Unear polarization, im- 
plying significant asphericity for its continuum-scattering envi- 
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ronment (Leonard et al. 2000). Prior to removal of the inter- 
stellar polarization, the polarization spectrum is characterized 
by a flat continuum (at p « 2%) with distinct changes in polar- 
ization associated with both the broad (symmetric, half width 
near zero intensity > 10,000 km s~') and narrow (unresolved, 
full width at half maximum < 300 km s"') line emission seen 
in the total flux spectrum. When analyzed in terms of a po- 
larized continuum with unpolarized broad-line recombination 
emission, an intrinsic continuum polarization of 3% re- 
sults, suggesting a global asphericity of > 45% from the oblate, 
electron-scattering dominated models of Hoflich (1991). The 
smooth, blue continuum evident at early times is inconsistent 
with a reddened, single-temperature blackbody, instead having 
a color temperature that increases with decreasing wavelength. 
Broad emission-line profiles with distinct blue and red peaks 
are seen in the total flux spectra at later times, suggesting a disk- 
hke or ring-like morphology for the dense (jie ~ 10^ cm"^) cir- 
cumsteUar medium, generically similar to what is seen directiy 
in SN 1987A, although much denser and closer to the progen- 
itor in SN 1998S. Such a disk/ring-like circumstellar medium 
may have formed from a merging of a binary-companion star 
that ejects the common-envelope material in the direction of 
the orbital plane of the binary system (Nomoto, Iwamoto, & 
Suzuki 1995). 

5. DISCUSSION 

5.1. SN1999em 

To interpret the X-ray and radio data on SN 1999em, we 
use the circumstellar interaction model proposed by Chevalier 
(1982), and elaborated by ChevaUer & Fransson (1994) and 



Fransson, Lundqvist, & ChevaUer (1996; hereafter, FLC96). 
The outer supernova ejecta are taken to be freely expanding, 
with a power law density profile 

p,n=At-\r/tr, (1) 

where A and n are constants. SN 1999em was a Type 11-P SN, 
so we expect that the star exploded as a red supergiant with most 
of its hydrogen envelope. The value of n for such a star is usu- 
ally taken to be in the range 7-12 (Chevalier 1982), although 
Matzner & McKee (1999) recently found that the explosion of 
a red supergiant leads to n = 1 1 .9 at high velocity. In the case of 
SN 1999em, we have some additional information from mod- 
eling of the optical spectrum. Baron et al. (2000) were able 
to model spectra on 1999 October 29 and November 4/5 with 
n = 7, 8, or 9. The modeling of the data for the earlier time sug- 
gests a relatively flat density profile (n = 7). The spectra on 1999 
October 29 also show evidence for a secondary absorption fea- 
ture at 20,000 km s~\ implying the presence of denser material 
in the cool ejecta. One possible origin for high velocity, dense 
ejecta is the shell that can form as a result of the diffusion of 
radiation at the time of shock breakout (Chevalier 1976); if this 
is the origin, it represents the highest velocity material in the 
ejecta. We assume that the outer density profile can be approxi- 
mated as a power law out to this high velocity and take n = 7 or 
9. The value of A is equivalent to specifying the density in the 
ejecta at some particular velocity and age. The model by Baron 
et al. (2000) for 1999 October 29 has p.,„ = 0.4 x lO'^" g cm'^ at 
10,000 km s^' at an age of 1 year. This value is in approximate 
accord with the density found in models with a broken power 
law density profile, mass of 10 M©, and energy of 1 x 10^' ergs 
(Chevalier & Fransson 1994). 
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The density in the wind is p„ = M /{4Trr^v„), where M is 
the mass-loss rate and v„ is the wind velocity. We take M-g = 
M/10~^ Mq yr~^ and v^i = v^/lO km s~^ as reference values. 
The interaction of the supernova with a freely expanding wind 
leads to a shocked shell with radius R oc t^-^(n = 7) or f^/^(n = 9). 
This shell is the source of the X-ray emission. The low observed 
X-ray luminosity suggests that the shocked gas is not radiative; 
this can be checked for consistency after a model is produced. 
Under these conditions, the X-ray luminosity is expected to be 
dominated by emission from the reverse shock region, which is 
relatively cool, as observed. The luminosity of the shell can be 
estimated from eq. (3.10) of FLC96 

L,^ = 2.0xl03^C(n-3)(n-4)V^V''-ii«/^' [^)' V,-' (l^f 

(2) 



the n = 7 case, which justifies our use of the adiabatic expression 
for the X-ray luminosity. The cooling time grows more rapidly 
than the age, so we expect non-radiative evolution through the 
time of our observations. If cooling were important, a dense 
shell would form that could absorb X-ray emission from the re- 
verse shock wave, as appeared to occur in the initial evolution 
of SN 1993J (FLC96; Suzuki & Nomoto 1995). Reradiation of 
the absorbed emission can give broad optical emission Unes, as 
observed in SN 1993 J at ages > 1 year. Such lines were not 
observed in SN 1999em, which is consistent with non-radiative 
evolution. 

From eq. (2.2) of FLC96 with n = 9, we find that the maxi- 
im velocity for n = 9 is 



ergs s 



where Q is 0.86 for solar abundances and 0.60 for n(He) / n(H) = 
1, = r/10^ K, and V4 is the peak velocity in units of 
10"* km s~'. Baron et al. (2000) find possible evidence for en- 
hanced He in SN 1999em, but this is uncertain and is probably 
not expected in a Type II-P SN. We take C, = 0.86; the uncer- 
tainty is small. On 1999 November 1, which we take to be 
day 4, the observed luminosity (0.4-8 keV) is 2 x 10^^ ergs s"^ 
and the temperature is 5.0 keV. With V4 = 1.5, we find that 
M-g/vwi « 1 (m = 9) or 2 (n = 7), quite similar to the value of 
M-g/vvvi 2 obtained above from the radio observations. With 
this mass-loss rate, the cooUng time for the gas, as deduced 
from eq. (3.7) of FLC96, is 31 days for n = 9 and is longer for 



15M 



1/7 



-1/7 



(3) 



where po specifies an ejecta density at a particular time to 
and velocity Vg, as described above. At / = 4 days with po = 
0.4 X 10"^° g cm"3 at Vo = 10,000 km s~> and to = 1 year and 
M-(,/vn.i = 1, we have V = 13,000 km s"'. For the n = 7 case 
with M-e/vy^i = 2, we have V = 14,000 km s"' and V oc 
These velocities are lower than the highest velocities deduced 
by Baron et al. (2000) on 1999 October 29 (day 1). However, 
there is some uncertainty in the high velocity, and there may 
have been rapid evolution at early times. 

The shock velocity determines the temperature of the reverse 
shock region. From eqs. (3.1) and (3.2) of FLC96, we have 
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Jrev = 2AVi keV for n = 9, so that with y4 = 1 .3, Jrev = 4. 1 keV. 
For n = 7, we have Trey = 4.8^4^ keV so Trev = 9.4 keV for 
V4 = 1.4. The observed temperature on 1999 November 1 falls 
between these two cases. The predicted evolution of the tem- 
perature is r oc y2 oc f-2/(«-2) oc („ = 9) and oc f""'* (n = 7). 
Table 3 compares this predicted evolution with that observed. 
The observations show a clear cooling, as expected. An addi- 
tional factor could be the lack of electron-ion equilibration at 
the reverse shock front; if it is not achieved, the electrons are 
cooler than the above estimate. From eq. (3.6) of FLC96, the 
expected conditions are such that the equihbration timescale is 
comparable to the age throughout the evolution, so that no firm 
conclusion can be drawn. The fact that cooling of the emission 
is clearly observed argues for a relatively flat density profile for 
the ejecta, as also found by Baron et al. (2000) from optical 
spectroscopy. 

An estimate of the expected evolution of the total X-ray lu- 
minosity, Lx, is given by Chevalier & Fransson (1994). When 
free-free emission dominates (r > 4 x 10^ K), Lx oc t~^, and 
when fines dominate (10^ < T < 4 x lO'^ K), cx r"-^^ (n = 9) 
or cx t~^-^^ (n = 7). The observed evolution is less steep than 
f"' and is closer to the expectation for line emission. The ob- 
served temperatures indicate that a transition occurred. Again, 
the observations are in approximate accord with expectations. 

One aspect of the X-ray observations that is not compatible 
with the simplest wind interaction model is the increase of the 
soft X-ray flux on 2000 February 2 (day 102, see Fig. 2). In 
our model, the emission is primarily from the reverse shocked 
ejecta, so the increase could be due to a region of increased den- 
sity in the ejecta. At this age, the velocity of the ejecta entering 
the reverse shock front is ^ 8,000 km s"'. Alternatively, the 
emission could be from an inhomogeneity in the circumstellar 
wind. The smooth wind is shocked to a temperature ^ 10^ K, 
so that a density contrast of ~ 100 is needed for the inhomo- 
geneity. The inhomogeneity must be large to affect the total 
luminosity, but a substantial change in stellar mass-loss charac- 
teristics at only 10'* cm from the star is unlikely. The fact that 
the change is transitory also argues against a change in mass- 
loss properties at this point. 

Although SN 1999em was detected as a weak radio source, 
the radio light curves are not well defined because of the very 
low radio flux density. The early observations were not of suf- 
ficient sensitivity to detect the source so the time of light curve 
peak and the early absorption process are also not well defined. 
The 8.44 GHz observations are the most useful in this regard. 
At late times, the observed flux evolution of other radio su- 
pemovae is oc or steeper (Weiler et al. 1986), so that the 
8.44 GHz upper limit on day 7 probably indicates that absorp- 
tion was significant at that time. With this value for the time 
of peak flux at 8.44 GHz (actually a lower limit), we can place 
SN 1999em on a diagram of peak radio luminosity vs. time 
of peak (Fig. 7). The diagonal fines show the velocity of the 
shocked shell if synchrotron self-absorption were the important 
process for the early absorption (Chevalier 1998). The posi- 
tion of SN 1999em gives about 10,000 km s"\ or less, as the 
shell velocity. Thus, it is just possible that synchrotron self- 
absorption was important at early times but more likely that an- 
other process, probably free-free absorption by the unshocked 
wind, was dominant. From eq. (2.3) of FLC96, the 8.44 GHz 
free- free optical depth on day 7 with V4 and M-^/v^i as de- 
duced above is r = 0.2T^ ^^'^ (n = 9) or 0.5T^ (n = 7), where 
75 is the temperature of the unshocked wind in units of 10^ K. 



The n = 7 case is somewhat more consistent with the observa- 
tions because of the higher mass-loss rate deduced for this case. 

The free-free absorption of the radio emission depends on 
the circumstellar temperature. Eastman et al. (1994) find from 
a simulation of a Type 11-P explosion an effective temperature at 
outbreak of Jeff = 1 -8 x 10^ K and luminosity 9 x 10'*'' ergs s"' . 
This temperature is similar to that used in one model by 
Lundqvist & Fransson (1991) for SN 1987A (their Table 1). 
Their model with M-^/vki = 0.4 gives a wind temperature of 
1 .2 x 10^ K for Teff = 2 X 10^ K. Because of the low wind den- 
sity, this temperature is not sensitive to the mass-loss rate and 
is mainly set by the outburst radiation temperature. An approx- 
imate temperature in the wind of SN 1999em should therefore 
be Ts « 1.2. The temperature at ^^30 days should not be very 
different because of the low wind density. X-ray heating by the 
shock has not been included here, but is probably not important 
because of the low X-ray luminosity. 

Overall, the X-ray and radio observations of SN 1999em 
are consistent with the non-radiative interaction of supernova 
ejecta with a power law profile (n ~ 7) interacting with a pre- 
supemova wind with M-g /v^i « 2. We estimate the accuracy of 
this mass-loss density to be a factor of 2. The X-ray emission 
is primarily from the reversed shocked ejecta, and the model 
can account for the luminosity, temperature, and temperature 
evolution of the emission. It does not account for the large soft 
X-ray flux observed on day 102. The radio observations are 
compatible with this general picture and suggest a comparable 
mass-loss rate. 

5.2. SN1998S 

The high X-ray luminosity of SN 1998S places it in the 
same class as the X-ray bright Type Iln SNe 1978K, 1986J, 
1988Z, and 1995N. Fox et al. (2000) discussed X-ray obser- 
vations of SN 1995N and summarized the results on the other 
SNe. SN 1995N had an X-ray luminosity of 1 .5 x lO'*' ergs s"' 
at an age of 2.0 years. SN 1998S is at the low end of the X-ray 
luminosities spanned by these sources 

The optical emission from SN 1998S shows both narrow and 
broad fines that can be attributed to circumstellar interaction 
(Fassia et al. 2000; Leonard et al. 2000; Chugai 2001). The 
narrow lines are presumably from dense circumstellar gas that 
is heated and ionized by the supernova radiation. Although 
high velocities (V4 « 1) are observed in the early evolution of 
SN 1998S, by day 494 the maximum observed velocities were 
V4 = 0.5-0.6 (Leonard et al. 2000). These velocities are much 
higher than the estimated velocity of gas at the photosphere, 
indicating that the emission is from circumstellar interaction. 

The interpretation of the X-ray emission from Type Iln su- 
pemovae is not straightforward because of the possibiUty that 
the emission is from shocked circumstellar clumps (see, e.g., 
Weiler et al. 1990 and Chugai 1993 on SN 1986J and Van Dyk 
et al. 1993 on SN 1988Z). The radio results for SN 1998S in- 
dicate that the CSM is likely clumpy or filamentary (Lacey 
et al. 2002). However, if our finding of heavy element over- 
abundances is correct, the X-ray emission is expected to be 
from shocked SN ejecta. If the interaction can be approxi- 
mately described by smooth ejecta running into a smooth cir- 
cumstellar wind, the temperature behind the reverse shock front 
is Trev = 2.4y2 keV (n = 9), 5AV^ keV (n = 7), or 9.6Vi keV 
(n = 6). It is difficult to produce the high observed temperature 
(r « 10 keV) unless there is unseen high velocity gas and the 
ejecta density profile is relatively flat. The observations do not 
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Fig. 7. — The peak spectral radio luminosity vs the product of the time of the peak and the frequency of the measurement, based on Fig. 4 of Chevalier (1998) with 
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flux density peak; an electron spectrum N(E) oc E~^-^ is assumed. 



show the cooling that is expected as the shock interaction re- 
gion decelerates, but the range in time is not large enough for 
the effect to be clearly visible. 

The X-ray luminosity shows a clear decline with time, ap- 
proximately (X t"^-^. This is in reasonable agreement with 
what is expected for free-free emission from a non-radiative 
shock region (f"' dependence as discussed above), so we can 
again use eq. (2) to estimate the mass-loss rate. The result is 
M_4/vwi = (1 -2), where M-4 = M/IQ-'^ Mq yr"', if n = 6 and 
¥4 = 1 (these estimates are impUed by the high observed tem- 
perature of '^lO keV) on day 665, again comparable to that 
determined from radio modeling of M-4 /v„i ~ 2. A high mass- 
loss rate is required to produce the X-ray luminosity, but the 
cooling time is still longer than the age. The assumption of a 
non-radiative reverse shock region is justified. In this model, 
the X-ray luminosity of the forward shock region is compara- 
ble to that from the reverse shock. However, the higher temper- 
ature of the forward shock region implies that the reverse shock 
dominates the emission in the energy band that we observed. 

The radio luminosity of SN 1998S is low for a RSN with 
a late turn-on time (see Fig. 7). This implies that synchrotron 
self-absorption was not important for the early absorption. This 
accounts for the basic features of the X-ray emission. The flat 
supernova density profile implies that much of the supernova 
energy is in high velocity ejecta. A similar deduction was made 
by Chugai & Danziger (1994) for SN 1988Z. The heavy ele- 
ment overabundances require that these elements be mixed out 
to a high velocity in the ejecta, which has been shown to occur 
in Cas A (Hughes et al. 2000). In theoretical models of aspher- 
ical explosions, heavy-element-rich matter is naturally ejected 
at high velocities, which well explains the peculiar late spectral 
features of SN 1998bw (Maeda et al. 2000). The polarization 
and spectral line profiles observed in SN 1998S (Leonard et al. 
2000) show that the explosion was complex, i.e. that the spher- 



ically syimnetric model described here must be an oversimpli- 
fication. 

We turn now to the observed elemental abundance ratios in 
the summed Chandra spectrum (Table 5), which are enhanced 
by a factor of ^ 3 - 30 over solar. These heavy elements must 
have been synthesized in the core evolution and explosion. 

In core-collapse supernovae, a strong shock wave forms 
and propagates outward through the star (e.g., Shigeyama & 
Nomoto 1990). Behind the shock, materials are explosively 
burned into heavy elements and accelerated to ~5000 km s~^ . 
At the composition interface, the expanding materials are de- 
celerated by the collision with the outer layers, which induces 
the Rayleigh-Taylor instability (e.g., Ebisuzaki, Nomoto, & 
Nomoto 1988; Arnett et al. 1989). Heavy elements are then 
mixed into outer layers. Such a mixing induced at the Si/O, 
C/He, and He/H interfaces brings Fe-peak elements, as well as 
Si, Mg, Ne, O, and C, into the H-rich envelope. These heavy el- 
ements in many Rayleigh-Taylor fingers are rather well-mixed 
through the interaction between the fingers (e.g., Hachisu et 
al. 1992). We thus assume that, although the heavy elements 
mixed into the H-rich envelope consists of only a fraction of 
the core materials, the abundance ratios among the mixed heavy 
elements are the same as in the integrated core materials. The 
mixing associated with the Rayleigh-Taylor instabilities contin- 
ues until the shock wave reaches the surface and the expansion 
becomes homologous. 

When the ejecta colhdes with the circumstellar matter, for- 
ward and reverse shock waves form and the latter propagates 
back into the H-rich envelope. The mass of the reverse-shocked 
H-rich envelope 6Mh at the time of the observations depends 
on the strength of the circumstellar interaction and the density 
structure of the H-rich envelope (e.g., Chevaher & Fransson 
1994; Suzuki & Nomoto 1995). We denote by SMz the mass 
of heavy elements mixed into the reverse- shocked H-rich en- 



SNel999emandl998S 



11 



velope. The observed mass fraction of heavy elements is then 
6Mz/XSMu, where X denotes the H mass fraction. 

We briefly summarize here the three types of heavy elements 
ejected from core-coUapse SNe (e.g., Thielemann, Nomoto, & 
Hashimoto 1996): 

O, Ne, and Mg: These elements are synthesized mostly 

before collapse during core- and shell-C burning; some 
Ne is processed into Mg and O during Ne burning. In 
the explosion, these elements are partially produced by 
explosive C burning, Ne is partially burned into Mg 
and Si, and O and Mg are partially burned into Si. The 
masses of these elements are sensitive to the progenitor 
mass, M, increasing with M. However, the degree of 
dependence on M varies across elements. For example, 
the Ne/Mg ratio depends on the conditions of C and Ne 
burning, thus depending on M. 

Si, S, and Ar: These elements are produced by explosive O 
burning during the explosion. Their masses do not 
much depend on M since the pre-coUapse core structure 
is not very sensitive to M. 

Fe: The mass of ejected Fe is very sensitive to the mass 
cut that divides the compact remnant and the ejecta, 
or the amount of fall-back matter. Theoretically, it 
is not possible to make a good estimate of the Fe 
mass because of the uncertainties in the explosion 
mechanism. 

In our comparisons, we use Si for normalization because it 
is the least dependent on M. We do not use Fe for comparison 
because of the large theoretical uncertainties. We now compare 
the O/Si, Ne/Si, and Mg/Si ratios (relative to the solar abun- 
dance ratios) with the theoretical calculations of Nomoto et al. 
(1997). These models are summarized in Table 6. 

The observed O/Si value of 0-1.2 seems to indicate that the 
mass of the progenitor was <20 Mq. The Mg/Si value of 0- 
0.7 also points to a progenitor mass of <20 Mq. The observed 
Ne/Si ratio of 0.6-14 excludes progenitor masses of <15 Mq. 
According to these models, the mass of the proginitor must have 
been between 15 Mq and 20 Mq. 

We note that the supernova yields are determined by the he- 
lium core mass Mue of the progenitor rather than the zero-age 
main-sequence mass M. Thus the Mg/Si and Ne/Si results ac- 
tually provide the constraint Mne = 5 ± 1 Mq (Nomoto et al. 
1997). The M-Mhc relation depends on the mass loss. If we 



adopt the spiral-in binary scenario for the origin of Type lln SN 
(Nomoto, Iwamoto, & Suzuki 1995), a large amount of mass- 
loss in the early phase of the progenitor's expansion leads to a 
smaller Mue compared with the case of no mass loss. In this 
case, the progenitor of SN 1998S could be initiaUy more mas- 
sive. 

6. SUMMARY 

Observations at radio, optical, and X-ray wavelengths have 
allowed the estimation of the physical parameters of two quite 
different Type II SN examples. Type II-P SN 1999em and Type 
Iln SN 1998S. From these results, it is possible to study, over 
a broad wavelength range, the physical parameters of the ex- 
plosions and the structure and density of the CSM established 
by the pre-SN wind. Such results are of importance for estima- 
tion of the properties of the pre-SN star and its last stages of 
evolution before explosion (e.g., Nomoto, Iwamoto, & Suzuki 
1995). In addition, the Chandra data have revealed the presence 
of many heavy elements in the spectrum of SN 1998S. We have 
shown the kind of analysis that can be done by comparing the 
observed X-ray determined abundance ratios to the model pre- 
dictions, and we were able to show that the mass of the progen- 
itor was ~ 15-20 Mq. We do caution, however, that our model 
assumptions on mixing are approximate and that there are un- 
certainties in both the interpretation of the X-ray data and the 
model calculations. Because of the great diversity of Type 11 
(and related Type Ib/c) SNe, only the multiwavelength study of 
additional examples can yield a comprehensive understanding 
of the last stages of massive star evolution. 
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Table 1 



Chandra Observations of SN 1999em. 



Date 


Dav 
( I'l'om Ref.) 


Length 
(ksec) 


Length After Filtering 

(ksec) 


Count Rate 
(cts ksec"') 


/x[2-10keV]/ 

J XL ' J' 

/J0.4-2keV] 


1999 Nov 1 


4 


29.0 


21.6 


3.70 


2.1 ±0.9 


1999 Nov 13 


16 


26.0 


10.6 


2.46 


0.8 ±0.4 


1999 Dec 16 


49 


35.3 


35.1 


1.02 


0.7±0.5 


2000 Feb 7 


102 


38.5 


23.9 


1.84 


0.2±0.1 


2000 Oct 30 


368 


26.4 


26.4 


0.46 


0.9 ±0.5 


2001 Mar 9 


495 


26.7 


24.6 


0.31 


0* 


2001 Jul 22 


630 


29.8 


25.1 


0.19 


0* 



*No flux detected above 2 keV. 



Table 2 

Chandra Observations of SN 1998S. 



Date 


Day 
(From Ref.) 


Length 
(ksec) 


Length After Filtering 
(ksec) 


Count Rate 
(cts ksec"') 


2000 Jan 10 


678 


18.9 


18.9 


34.4 


2000 Mar 7 


735 


23.4 


23.1 


30.3 


2000 Aug 1 


882 


19.8 


19.8 


25.8 


2001 Jan 14 


1048 


29.2 


29.2 


19.5 


2001 Oct 17 


1324 


28.7 


14.8 


12.4 



Table 3 



X-RAY Temperature Evolution of SN 1999em. 


Day 
(From Ref.) 


Obs. MEKAL kT 
(keV) 


Obs. Brems. kT 
(keV) 


Model kT {n = 1) 
(keV) 


Model kT (n = 9) 
(keV) 


4 

16 

49 

102 

368 


5.0!^:! 

2 5+2-7 

^•-'-0.9 

2.9+2-^ 
0.8!«i 


5 2+"-5 

9 a+5.9 
^•-'-1.1 
2 7+3.2 

^•-'-0.4 
1 c+6.5 
^••^-0.8 


9.2 
5.4 
3.4 
2.5 
1.5 


4.1 
2.8 
2.0 
1.6 
1.1 



Note. — Listed uncertainties are 90% confidence intervals. 



Table 4 

VMEKAL Fits to the Spectra of SN 1998S. 



Day 
(From Ref.) 


kT 
(keV) 


Fe Abund. 
(w.r.t. solar) 


X^/d.o.f. 


Lx[2-10 keV]'' 
(ergs s"') 


678 
735 
882 
1048 


io.4!ii/ 

10.4!iV 


9 0+15.8 
2 8 
5 4+14.5 


82.7/96 
71.7/71 
65.8/75 
57.0/59 


9.3 X 10^^ 
8.8 X 10^' 
6.1 X 10^' 
5.3 X 10^9 



Note. — Listed uncertainties are 90% confidence intervals. 
^forz = 0.003 and7/o = 50 
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Table 5 

Elemental Abundances from VMEKAL Fit to Summed Spectrum of SN 1998S. 



Element 


Best-fit Abund. 


90% Confidence Interval 


O 


0.7 


0-2.9 


Ne 


15 


8.7-35 


Mg 





0-1.6 


Al 


33 


0-135 


Si 


5.7 


2.4-15 


S 


8.7 


3.0-24 


Ar 


18 


3.4-52 


Ca 


0.8 


0-13 


Fe 


3.0 


1.8-6.8 


C 


1.0 


0-200 


N 





0-18 


Na 





0-36 


Ni 





0-14 



Note. — These abundances are actually the ratio of a 

given element to H normalized to the similar solar ratio. 
For example, "Ne" is actually [(Ne/H)98sJ/[(Ne/H)0]. 



Table 6 

Elemental Abundance Ratios Relative to Solar for SN 1998S and Theoretical Models. 



Abund. ratio 


SN 1998S 


13 Mq 


15 Mq 


18 Mq 


20 Mg 


25 Mq 


40 Mq 


Ne/Si 


0.6-14 


0.14 


0.17 


0.86 


1.1 


2.2 


0.56 


Mg/Si 


0-0.7 


0.18 


0.49 


0.58 


1.09 


2.1 


1.1 


O/Si 


0-1.2 


0.20 


0.43 


0.80 


1.4 


2.4 


1.6 



